Hole Dynamics in Two-Dimensional Antiferromagnetic Mott Insulators 
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The dispersion relation of a doped hole in the half-filled 2D Hubbard model is shown to follow a 
\k\ 4 law around the (0, ±7r) and (±7T, 0) points in the Brillouin zone. Upon addition of pair-hopping 
processes this dispersion relation is unstable towards a \k\ 2 law. The above follows from T — 
Quantum Monte calculations of the single particle spectral function A(k, u) on 16 x 16 lattices. We 
discuss finite dopings and argue that the added term restores coherence to charge dynamics and 
drives the system towards a d x 2_ y 2 superconductor. 



PACS numbers: 71.27.+a, 71.30. +h, 71.10.+X 

The excitation spectrum of a band insulator is gapped 
in both the spin and charge degrees of freedom. In this 
case the single particle Green function will essentially 
take a free particle form since there are no low lying ex- 
citations on which the added particle can scatter. The 
Mott insulator with long-range antiferromagnetic order 
has a charge gap but no spin gap. This leads to non- 
trivial hole dynamics which will depend on the dynamical 
spin susceptibility as well as on the coupling between spin 
and charge degrees of freedom. Hole dynamics in mag- 
netic insulators have been addressed in the pioneering 
work of Brinkman and Rice [jjj . Progress in photoemis- 
sion spectroscopy has provided us with an experimental 
realization of this problem In this Letter, we show nu- 
merically that hole dynamics in a two-dimensional Mott 
insulator with long-range magnetic order may behave in 
radically different ways. The dispersion relation of a 
hole doped into the half-filled Hubbard model is very flat 
around the (0, ±7r) and (±7r, 0) points in the Brillouin 
zone. We will give numerical evidence that it follows a 
| A; | 4 law. We argue that this flatness and resultant singu- 
lar momentum dependence of charge excitations provides 
us with a microscopic basis for understanding incoher- 
ent charge dynamics and unusual character of the metal- 
insulator transition. Upon addition of a term with matrix 
element W which describes pair-hopping processes, this 
flat dispersion relation is unstable towards a \k\ 2 law. Al- 
though the added term for small values of W does not 
alter the insulating state itself it restores coherence to 
charge dynamics in the vicinity of the Mott transition. 

The model we consider reads: 
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component of spin a on lattice site i and 

1/2. We consider a square lattice of linear size L and 
impose periodic boundary conditions in both lattice di- 
rections. The t-U-W model has been introduced and 
studied in Ref. ||-|§ . The results show that at half- filling 
and finite values of U/t, W triggers a quantum transition 
between the Mott insulator and a d x 2_ y 2 superconductor. 

To understand the influence of the W-term on hole 
dynamics it is convenient to consider the strong cou- 
pling limit where double occupancy is prohibited. In this 
limit, the W-teim describes pair-hopping processes in the 

spin singlet, -W^s/' ( A ?,A/' + A ?,5 A ?/') and Spin 
triplet W£. (tL T->* + T ? * tJ, ) chan- 
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creates an electron with z- 



TJ- creates a spin triplet on the bond i, i + 5 with 
^-component of spin m. In the strong coupling limit, 
the states TJ - and Al - form a complete basis for a 

i,S,m i,8 

pair of electrons on the bond i, i + S. Thus, any spin 
configuration surrounding a hole may be written as a 
superposition of triplets and singlets on a set of bonds 
covering the lattice. One may now see how hole dynam- 
ics are effected by the W-term: the motion of a pair of 
electrons is nothing but the hopping of a hole within the 
same sublattice accompanied by a local rearrangement of 
spins. Due to zero point quantum fluctuations of the spin 
background, the resultant state is not orthogonal to the 
ground state. Since the state in which the hole is doped 
has antiferromagnetic order, we expect the singlet pair- 
hopping processes to dominate the low-energy physics. 
Those processes are similar to the three site term ob- 
tained in a strong coupling expansion of the Hubbard 
model 0. The influence of those terms on hole dynam- 
ics has been addressed in Ref. [||. 

At half-band filling auxiliary field quantum Monte 
Carlo (QMC) simulations provide an efficient tool for 
the study of the above model. Partly due to the pres- 
ence of particle-hole symmetry, the infamous sign prob- 
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lem is avoidable. Our T = data are produced with the 
projector QMC (PQMC) algorithm @0 supplemented 
with a numerically stable method to compute imaginary 
time displaced correlation functions JO]. The imaginary 
time data is analytically continued to the real axis with 
the Maximum- Entropy method |l^|l3|]. We have used 
a flat default model and taken into account correlations 
in imaginary time data with the use of the covariance 
matrix. 

We start by considering the half-filled Hubbard model 
at U/t = 4 and T — 0. For this parameter set and 
after extrapolation to the thermodynamic limit numeri- 
cal simulations lead to the conclusion that the ground 
state is an insulator with single particle gap A qp = 
0.67 ± 0.015 (TJ] and long-range antiferromagnetic or- 
der: to = liniL^oo \J 3(m z ((L/2, L/2))m z (6)) = 0.39(5) 
Here, m z (i) = - n^. Fig. |l| plots 

A(k, ui) = lmG(k,ui), where G(k,u>) corresponds to the 
single-particle Green function. Due to particle-hole sym- 
metry, A(k, oj) = A(k + Q, —uj) with Q — (tt, tt) in units 
of the lattice constant. The sum rule: J_ oa A(k,oj)doj 
= tt (ct c^r a ) is satisfied. We are primarily in- 
terested in the single-hole dispersion relation E(k) as 
defined by the peak position in A(k,uj). Around the 
(0, tt) and three equivalent points in the Brillouin zone, 
E(k) shows a very flat structure. Along the (0, 0) to 
(0, 7r) direction E(k) is compatible with a \k\ A law over 
an acceptable range (see Fig. |). The same conclu- 
sion is reached when considering the (tt, tt) to (0, tt) di- 
rection. Such a flat dispersion relation has been ob- 
served numerically in t-J model calculations ]l7p . The 
overall flatness appears consistently with the flat disper- 
sion observed in underdoped cuprates . The fol- 
lowing points are of interest, i) The energy difference 
AE = E(k = (tt/2, tt/2)) - E(k = (0,tt)) is not un- 
ambiguously distinguishable from zero within our reso- 
lution. AE = 0.045i, 0.06i and 0.015t for the L = 8, 12 
and L = 16 lattices respectively. The uncertainty of our 
data, is in the same ball-park as the above quoted val- 
ues. We note that the flatness of the dispersion relation 
around the (0, tt) point should lead to a broader lineshape 
at (0, tt) than at (7r/2, tt/2) [^0| thus leading to some am- 
biguity in the definition of E(k). In fact, defining E(k) 
as the leading edge rather than the peak position yields a 
negative value for AE for the L = 16 lattice. Lineshapes 
are notoriously hard to compute with the Maximum- 
Entropy method and further large scale calculations are 
required to confirm this point, ii) A shadow-band due to 
the presence of long-range magnetic order is seen along 
the (tt/2, tt/2) to (tt, tt) to (tt, 0) direction. Around the 
(0, 0) point a two peak feature is seen in the data. A 



ters constant, U/t = 4, T = 0, (n) = 1. At this point 
in parameter space, the ground state remains a Mott in- 
sulator with to = 0.24(1) and A qp = 0.54(6) g. Fig. 
| plots A(k, u) again on a L — 16 lattice. The follow- 
ing points are of importance, i) Upon inspection one 
notices that the bandwidth is substantially enhanced by 
the inclusion of the VF-term. This high energy phenom- 
ena may be captured at a mean-field level by the Ansatz: 
(n ?jT ) + (n ?4 ) = 1, (n ?>T ) - (n ?4 ) = (-m HF ) l *+ l y and 
(Kfi = K. Self-consistency yields the single particle 
gap Ahf = Vmnpj2 = 0.43t and a band ranging from 
— 11.9£ to 11.9t. The bandwidth agrees well with the 
Monte-Carlo data. This approximation underestimates 
the single particle gap thus showing that is does not cap- 
ture the low-energy physics contained in the W-tevm. 
ii) Around the (0, tt) point, E(k) is not as flat as for 
the Hubbard model. As shown in Fig. 0, the data is 
compatible with a quadratic fit. To be more precise, 
as W/t is enhanced the domain in fc-space around the 
(0, ±7r), (±7t, 0) points which is compatible with a \k\ 4 fit 
is suppressed, iii) As in Fig |l| the shadow-band feature 
is present, iv) The energy difference in the dispersion 
relation between the (tt/2, tt/2) and (0, tt) points is not 
distinguishable from zero. |22j] 

Next we study finite dopings. Here, we are confronted 
to a sign problem so that the CPU time required to 
achieve a given precision scales exponentially with inverse 
temperature (f3) and lattice size. The presented data 
is produced with the finite temperature QMC method 
p3| , ^6[ . We first consider the vertex contribution to pair- 
ing correlations in the d-wave and extended s-wave chan- 
nels. This quantity is defined by: 

Pt s (r) = (At s (r)A d , s (0)) - ]T f d , s (8) f d , 3 (6') (3) 



similar feature at larger couplings is seen in Ref. |21|. 
We now set W/t = 0.15 and keep the other parame- 



Where A L( f ) = T, a J fd,s($W C r,a C l + l_ a > ^ = 1 

fd(5) = 1(— 1) for 5 = ±(? x (±a^). Per definition, in the 
absence of interaction P^ s vanishes (i.e. in the above 
mentioned mean field approximation which takes into 
account band width effects, P^ s vanishes.) Due to the 
strong variation of P$ s (0) (24j] as a function of doping, we 

consider the quantity: P% s (r) = ^ s (^/^s(") which 
measures the decay rate. This quantity is plotted ver- 
sus doping in Fig. |]. At the largest distance on our 
8x8 lattice, the tU-term substantially enhances the d- 
wave signal. We note that the same conclusion is reached 
when considering P djS = (A d s (r)A d}S (0)). 

The uniform spin susceptibility is plotted versus tem- 
perature in Fig. ||(a) a t ( n ) = 0.78, U/t = 4, for var- 
ious values of W/t. Starting from the noninteracting 
case (U=W=0) and turning on the Coulomb repulsion 
to U/t = 4 enhances the spin susceptibility. At high tem- 
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peratures, this enhancement may be understood within 
the random phase approximation. As W/t grows, there 
is a suppression of the spin susceptibility. There are two 
effects which cause this suppression, i) The enhancement 
of the bandwidth as a function of W/t. ii) The growth 
of d-w&ve pair correlations. Disentangling the contribu- 
tion of those two effects at different energy scales is non- 
trivial. The spin structure factor S(q) = (m z (q)m z (—q)) 
is plotted in Fig. || (b). The reduction of S(q) at q = 
as W/t is enhanced may be traced back to our discus- 
sion of the spin susceptibility since the latter quantity is 
given by (3S(q = 0). In the vicinity of q = S(q) 
shows a somewhat sharper feature at W/t = 0.15 than 
at W/t = 0, thus showing that the magnetic length scale 
is enhanced by the W-teriri. 

In summary, we considered the t-U-W model at U/t = 
4 and for two different choices of W: W/t = 0, W/t — 
0.15. In both cases, the ground state at half-filling is a 
Mott insulator with long-range antiferromagnetic order. 
However, the nature of hole dynamics in this antiferro- 
magnetic background is strongly dependent on the choice 
of W/t. In the case of the Hubbard model we concluded 
that the hole dispersion relation is consistent with a \k\ 4 
law around the (0, ±7r) and (±tt, 0) points in the Bril- 
louin zone. This flat dispersion relation is compatible 
with the picture of incoherent charge dynamics and intro- 
duces a singular momentum dependence in the electron 
self-energy. Recently, the authors of Ref. |2^,^(| have 
computed the Drude weight on t-J clusters and found 
results consistent D ~ <5 2 . The sum rule for the optical 
conductivity is proportional to the doping 6 in the case 
of the t-J model. Thus, close to the metal insulator tran- 
sition, only a vanishingly small portion of the weight in 
the optical conductivity will be contained in the coherent 
Drude response. Introducing the TF-term alters this situ- 
ation. On one hand we have shown here that the disper- 
sion relation around the (0, ±7r) and (±7r, 0) points follow 
a \k\ 2 law. On the other hand, the Drude weight satisfies 
D ~ S for the t-J-W model jHH]. Thus, the IF-term 
restores coherence to charge dynamics in the vicinity of 
metal-insulator transition. However, the short-range an- 
tiferromagnetic spin correlations at finite doping remain 
robust upon switching on W. In the doped phase, d-wave 
pairing correlations functions are substantially enhanced 
by the inclusion of W thus lending support to the occur- 
rence of a superconducting state. In terms of quantum 
phase transitions the inclusion of the VF-term alters the 
dynamical exponent from z = 4 to z — 2 p^,p|, ^5| , p6[ . 
More generally, the W-tenm exploits one of the many 
potential instabilities of the incoherent metallic state re- 
alized in the Hubbard model in the vicinity of the metal- 
insulator transition. As a more realistic model for high- 
T c cuprates, smaller values of W are to be considered 
so as to study the interplay between the pairing energy 
scale and the larger energy scale associated with the flat 



bands. 

The numerical simulations were performed at the Su- 
percomputer Center of the Institute for Solid State 
Physics, University of Tokyo, the HLRS-Stuttgart and 
HLRZ-Jiilich. We thank M. Brunner, A. Muramatsu, H. 
Tsunetsugu and M. Zacher for conversations. 



A(k,u) 



FIG. 1 

at U/t = 




A(k, w) at T = for the half-filled Hubbard model 
4 on a 16 x 16 lattice. The considered path in the 
Brillouin zone is listed on the left hand side of the figure. We 
have normalized the raw data by the factor listed on the right 
hand side of the figure. This normalization sets the peak value 
of A(k, uf) to unity for all considered k vectors. 
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A( k,uj), (n) = 1, U/t = 4, W/t = 0.15, T = 0, L 




C//i = 4, £ = 8, /? = (L/2.L/2) 



2. Same as Fig. M but for W/t = 0.15. 



(n) = 1,L = 16,T = 0, U/t = 4 




(0,3tt/8) (0,tt) (0,37r/8) 

FIG. 3. Dispersion relation, E(k), as defined by the peak 
position in A(k,u). The solid lines correspond to a least 
square fit to the form a + b(n — k y ) 2 (a + b(-K — k y ) 4 ) for 
W/t = (W/t = 0.15). For the fit we consider the k y range 
[5n/8,n] for both choices of W/t. 
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FIG. 4. Equal-time pair field correlations versus filling in 
the d-wave and extended s-wave channels. s (R) = PJ" s (R) / 
P$ g (Q). The same scale is used for both (a) W/t = 0.15 and 
(bj W/t = 0. 



U/t = 4, {n) = 0.78, L - 
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(0,0) (0,7r) (it, it) 
T/t 3 

FIG. 5. (a) Spin susceptibility versus temperature for var- 
ious values of W/t and (n) = 0.78. (b) Spin structure factor 
at (3t = 6, (n) = 0.78 for the Hubbard and t-U-W mod- 
els. In both (a) and (b) the solid line corresponds to the 
non-interacting (U = W = 0) case. 
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